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Abstract 

We investigate the effects of a baryon-antibaryon asymmetry on the spec- 
trum of dileptons radiating from a quark gluon plasma. We demonstrate the 
existence of a new set of processes in this regime. The dilepton production 
rate from the corresponding diagrams is shown to be as important as that 
obtained from the usual quark-antiquark annihilation. 
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I. INTRODUCTION 



Experiments are now underway at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven to study nuclear collisions at very high energies. The hope is to produce a 
plasma of deconfined quarks and gluons. This plasma is, however, rather ephemeral and 
soon hadronizes into a cornucopia of mesons and baryons. One thus needs an indirect 
means of deducing as to whether or not the plasma was produced in the history of a given 
collision. Various experimental signatures have been proposed to this effect: J/^' suppres- 
sion 0], strangeness enhancement 0, dilepton spectra 0-^1 etc. In this paper we calculate 
a new contribution to the spectrum of dileptons ( i.e., e~^e^ , ^^^~) emanating from a quark 
gluon plasma. 

Early calculations of the dilepton radiation in the deconfined sector were concerned with 
the process qq —>■ e^e~ A recent calculation has estimated the effects of chemical 

non-equilibrium and of a large gluon excess on dilepton spectra There, a fugacity 
was introduced to account for chemical non-equilibrium. The function of this fugacity is 
essentially to change the gluon and quark numbers from their equilibrium values. Possible 
sources of dileptons, such as q + q annihilation, q + g Compton scattering, and g + g fusion 
had been investigated and at chemical and thermal equilibrium the spectrum was found to 
be dominated by qq e~^e~ , followed by qg — >■ qge~^e~ which is an order of magnitude lower, 
followed by gg qqe^e~ which is lower than the first process by 3 orders of magnitude 0. 

The aim of our work is to propose that, when there is an asymmetry in the populations 
of quarks and antiquarks (i.e., a finite baryon chemical potential) a new set of diagrams 
actually arise. Using these we calculate a new contribution to the 3-loop photon self-energy. 
The various cuts of this self-energy contain higher loop contributions to the usual processes 
of qq e~^e~, qg qe~^e~, qq — > qqe~^e~, and an entirely new process: gg e^e~ . We 
calculate the contribution of this new channel to the differential production rate of back-to- 
back dileptons. It is finally shown that within reasonable values of parameters this process 
may become larger than the differential rate from the standard tree level qq — >■ e^e~ . 

Imagine a scenario where the plasma is not just heated vacuum, but actually displays 
an asymmetry between quarks and antiquarks. This asymmetry would eventually manifest 
itself as an asymmetry between the baryon antibaryon populations in the final state. In 
plasma calculations this asymmetry may be achieved by the introduction of a quark chemical 
potential /ig. For the sake of simplicity we will assume here that is the same for u, and 
s quarks (we are assuming the plasma to contain three massless flavours). We also assume 
that the chemical potential for gluons is zero. The rest of the paper is organised as follows: 
section II discusses a class of diagrams which are non-existent at zero temperature, and 
also at finite temperature and zero density. These become finite at finite density. Section 
III focuses on a specific channel which will become a source of dileptons. In section IV we 
derive a new contribution to the photon self-energy at three loops and discuss its various 
cuts. Finally, in section V we calculate the production rate of back to back dileptons. All our 
calculations are done in the imaginary time formalism of equilibrium thermal field theory 
(our notation is described in Appendix A). We have assumed the presence of only three 
massless flavours of quarks. 
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II. NEW DIAGRAMS FROM BROKEN CHARGE CONJUGATION 

INVARIANCE 



At zero temperature, and at finite temperature and zero baryon density, diagrams in 
QED that contain a fermion loop with an odd number of photon vertices (e.g. Fig. |lD are 
cancelled by an equal and opposite contribution coming from the same diagram with fermion 
lines running in the opposite direction (Furry 's theorem This statement can also be 

generalized to QCD for processes with two gluons and an odd number of photon vertices. 




FIG. 1. Diagrams that are zero by Furry 's theorem and extensions thereof at finite tempera- 
ture. These become non-zero at finite density. 



However, at finite density or alternatively at finite fermion chemical potential, this can- 
cellation no longer occurs. As an illustration consider the diagrams of Fig. ^ for the case of 
two gluons and a photon attached to a quark loop (the analysis is the same even for QED 
i.e., for three photons connected to an electron loop). In order to obtain the full matrix 
element of a process containing the above as a sub-diagram one must coherently sum con- 
tributions from both diagrams which have fermion number running in opposite directions. 
The amplitude for r^'P''{= Tf'P'' + T^p^") are : 

{q + p-k)aqp{q + p)s 
{q + p- kYq^iq + pY ' 



and 



{q + k-p)o,qp{q~p)5 
{q + k — pYq^{q — pY 



(2.1) 



At finite temperature (T) and density (chemical potential ju), we have the zeroth com- 
ponent of the fermion momenta given by. 
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go = i{'2n + l)7iT + /i W n E I. 



(2.2) 



We assume that /i is the same for both flavours of quarks. Note that the extension of 
Furry's theorem to finite temperature does not hold at finite density: as, if we set n — *• 
—2n — 2 we note that go ~/*~ Qo and result 

T^P'i^i.T)^ -T'^P^'i^i.T). (2.3) 

Of course, If we now let the chemical potential go to zero (/i — * 0) we note that for the 
transformation n — >■ —2n — 2 we obtain go —go and thus T^^'^(0,T) — > — T'^''^(0, T). The 
analysis for fermion loops with larger number of vertices is essentially the same. 

A more physical perspective is obtained by noting that all these diagrams are are en- 
countered in the perturbative evaluation of Green's functions with an odd number of gauge 
field operators. At zero (finite) temperature, in the well defined case of QED we observe 
quantities like (0|A^iA^2---^/x2n+i|0) ( Tr[p(/i, /?)A^iA^2---^/x2ri+i] ) under the action of the 
charge conjugation operator C. In QED we know that CA^C^^ = —A^. In the case of the 
vacuum |0), we note that C|0) = |0), as the vacuum is uncharged. As a result 

(0|A^lA^2...AM2n+l|0) = {Q\C-^CA^,C-^CA^2.-A^2n+lC-^m 

= (0|A^iA^2...A^2n+l|0)(-l)'"+^ 

= -(0|A^lA^2-^M2n+l|0). (2.4) 

Hence, all such Green's functions are zero. The argument is the same for the case of finite 
temperature and zero density (or chemical potential). In the case of nonzero chemical 
potential, however, the system is charged, hence the eigenstates of the density operator p 
are not eigenstates of C i.e., C\n) ^ \n). The medium, being charged, manifestly breaks 
charge conjugation invariance and these Green's functions are thus finite. The appearance 
of processes that can be related to symmetry-breaking in a medium has been noted before. 



III. THE TWO-GLUON-PHOTON VERTEX 

Let us now focus our attention on the diagrams of Fig. ^ Such a process does not exist 
at zero temperature or even at finite temperature and zero density. At finite density this 
may lead to a new source of dilepton or photon production. In this section we obtain general 
expressions for this diagram. 

We sum the Matsubara frequencies using Pisarski's non-covariant method [|T^] (Our 
notation is however different from that in and is explained in Appendix A.). In this 



method, propagators may be expanded as follows 

1 ~ f-iP 

Q 

where 



1 = A(g) = {-^) r^dxl e-^^^'-^>' K'{\q\, -x"), (3.1) 
q-^ Jo 



S2 
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FIG. 2. The two gluon photon effective vertex as the sum of two diagrams with quark number 
running in opposite directions. 



where = 1 — n{E + fj,) , = —n{E — /x), and n is the Fermi-Dirac distribution function. 
Then, 

A{q + p — k) 



{q + p- k) 
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i-i) [ dx\ e^(9W-fc"-MK A(|g + p- (3.2) 



K{\q + p-klxl) = ^^^/,^(E^+^_,,.,;,)e-^(W.-iMK, 

where /+ = 1 — h{E — /i) , /_ = —h{E + /x). For the last propagator, we may write 
1 



A{q + p) 



{q+pf 

= {-%) dxl c^C'jW-mM A(|g + pl,x°) (3.3) 

JO 



with 



1 

In the above three equations Si, S2, S3 are sign factors which are either +1 or -1, and g° — *■ 
— ifi). We substitute the above three equations in Eq. (|2.1| ) and add the two diagrams. 
Now, one may simply perform the Matsubara sum to obtain a delta function over the times 
(see Appendix B). Using the delta function we may evaluate the remaining time integrations 
to give the full expression for this vertex as (see Appendix B) 

J (2vr)3 2 J 7 7 7 7 7 J ^0 + ,^e, - s,E,^, 
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+ SlEg+p-fc - S^Eq+p 

S2{h{Eg+p_k + sifi) - h{Eq+p_k - sin)) - si{n{Eq + ga/^) - n{Eq - S2fi)) 



(3.4) 



In the above equation Eg^k = ylo]'^ + I^P — 2|g||A;| cos6', where cos6' is the angle between q 
and k. While in the numerator 



Qs2 = (s2, sin ^ COS0, sin ^ sin 0, cos( 



and 



q sin 6 cos (p q sin 6 sin g cos 6 — k 



k\ 



Up to this point no approximations have been made. Completing the remaining angular 
integrations will lead to the expression for the effective vertex with arbitrary momenta k, p. 
This, however, turns out to be rather difficult. Even in the Hard Thermal Loop (HTL) 
approximation, evaluation of similar three point loops for arbitrary three momentum p is a 
difficult problem |T^. For the HTL approximation to be valid, the temperature should be 
high, and the momentum in the outer legs should be low. Predicted, initial temperatures 
at RHIC and LHC range from 300 to 800 MeV |13,14|. Dileptons from the plasma are 
predicted to become important in the intermediate mass regime (1-3 GeV) [^. In this region 
of parameter space the HTL approximation obviously does not apply . Hence we do not 
make this approximation in this work. Furthermore, in the interest of technical simplicity, 
we proceed to evaluate the above diagram in the limit of the photon three momentum p = 0, 
i.e., for dilepton pairs produced back-to-back. 

In this limit the expression for the effective vertex reduces to 



fipu 



(Pq & 



ah 



ssAnjEq^k, sijj) - siAn{Eq, g3/i) _ 



+ S2Eq - s^Eq 
S2Ah{Eq_k, sifi) - siAh{Eg, S2fi) 
k^-p^ + SiEq^k - S2Eg 



(3.5) 



In the above, the notation Ah{E, s/i) = h{E + s/i) — h{E — s/i)) is introduced. The An's 
quantify the effect of a finite chemical potential, as they represent the difference between 
the distribution functions of a quark and its antiquark. If — s> 0, all the An's go to zero 
and so does the entire expression. Noting that S2An{Eq, s^fi) = SsAn^Eg, Sa/i), we may add 
up the two terms in the integrand to give. 



•' l^^J ^ S1S2S3 

S2SiAn{Eg, ssfi) - S2S-iAn{Eg^k, Sifx) 

(fcO + S,Eg_k - S3Eg){kO - / + S,Eg_k ' ^2^, 



(3.6) 
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Note that the integrand of the above expression for the vertex has factorized into a function 
of the form /i(fco)/2(Po ~ ^o)- This result will become important in the eventual evaluation 
of the photon self-energy in the subsequent sections. 

We note that the sole dependence is in the numerator. This allows us to integrate 
out the (p dependence. The presence of the Lorentz indices a, f3, 7 indicate that Eq. ( p.6|) 
actually contains 64 separate terms. The Tr [7'^7'^7^7'^7'^7°] combines these into 64 other 
terms. Ignoring the trace part, the rest of the formula is found to contain only 20 nonzero 
terms (see Appendix C). Of these 20 terms we observe many to be interdependent. There 
are at most 8 independent terms (see Appendix C). 

To perform the 6 integration we shift q in the second term of the integrand of Eq. ( p.6|) 
to give 



^ J (27)3 2^ Tr[YTYl 7 7 



S1S2S3 



, . . S2Sl^n{Eq, 53/i) 



■ (3.7) 



The resulting 



integration now becomes rather simple. 



IV. THE PHOTON SELF-ENERGY AND ITS IMAGINARY PART 

We are now in a position to calculate the contribution made by the diagram of Fig. |^ to 
the dilepton spectrum emanating from a quark gluon plasma. To achieve this aim we choose 
to calculate the imaginary part of the photon self-energy as represented by the diagram of 
Fig. ^. In the previous section we wrote down expressions for T^P^{k —p,k]p). To write 
down the expression for the full self-energy we also need expressions for T^'"'(fc, k — p; —p). 
A simple analysis consisting of essentially reversing the direction of the internal momentum 
q leads us to the result that T'''"'(fc, k — p; —p) = T'^^'^{k — p, k; p). 

We may write down the full expression for the photon self-energy as 

1 r d^k 

K = ^Y.j j^,V,^{k)T^"'{k-p,k-p)V,^{k-p)Ti'>^{k,k-p;-p). (4.1) 

We perform this calculation in the Feynman gauge for the gluons, thus 

-D.pik) = ||^. (4.2) 

We calculate in the limit of photon three momentum p = 0. We now shift k —>■ k+p. Notice 
that, in the limit p = 0, this only implies shifting k^ by p^. This is followed by switching 
— > ~k^. Note that neither of these operations has any effect on the full nature of the 
integral as k^ is summed over all even discrete frequencies. With this we now obtain the full 
self-energy as 
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FIG. 3. The Full photon self-energy at three loop and the cut that is evaluated in this paper. 



'-'-p — y g '-'a/37 



d^kd^qdq'^ 1 



(2vr)^ 



E E 



X 



X 



q'-k 



s'^E,.) 



(4.3) 



Where the multiplicative factor C'^p^'^ is the result of the contraction of Lorentz indices 
appearing in the two traces of 6 7 matrices and the intervening factors of the metric. In the 
above expression, we have also assumed the presence of three massless flavours of quarks. 
For simplicity the chemical potential is assumed to be the same for all three flavours. In 
order to calculate the differential rate of back to back dileptons we need to evaluate the 
imaginary part of this diagram. This may be obtained in two ways: the conventional 



method involves converting the sum over discrete frequencies into a contour integral, 
followed by evaluation of the contour integral by summing over its residues and then looking 
for poles and branch cuts in the final expression in terms of p° by analytically continuing 
p° onto the real axis. Ten years ago Braaten, Pisarski, and Yuan presented an identity 
which elegantly achieved the end result of this procedure for a fermionic Matsubara 



frequency k^ . In the following we will refer to this formula as the BPY formula. It is given, 
generally, in terms of products of arbitrary functions of discrete imaginary frequencies as 
(for completeness, a brief proof of the BPY formula for a Bosonic frequency in terms of 



8 



the method of |TT| is presented in appendix D) 

/•+00 r+oo 

DiscT^/i(A;°)/2(p° - = 27rz(l - e^/^) / du / du' 

n{uj)n{uj')5{E - uo - uj')pi{uj)p2{uj'), (4.4) 

where pi(uj),p2(uj') are the spectral densities of /i(z),/2(z). Note that Eq. (^4.3|) is precisely 
in the form required by Eq. ( [4. 4]) , with 

1 



[pO-A;0-(siX-S3g)][A;o-(s'ix'-s!5g')][(P°-'^°)'-'^']' 
where, for brevity we have x{x') = Eq_k{Eqi^k) and g(g') = Eq{Eqi). From here, one may 
write down the spectral densities of the two functions by mere inspection: 

Pi{uj) = 6{uj - k)Res[fi{uj = k)] + 6{uj + k)Res[fi{uj = -k)] 

+ 5{uj - {siX - S2g))Res[/i(cj = SiX - S2q)] 

+ 6{uj - {s[x' - 4g'))Res[/i(cu = s[x' - s'^q')] 
P2{uj') = 6{uj' - k)Res[fi{uj' = k)] + 6{uj' + fc)Res[/i(cu' = -A;)] 

+ 6{uj' - {siX - S3q))Res[fi{uj' = SiX - Sag)] 

+ 6{uj' - {s[x' - s'sq'))Res[fi{uj' = s[x' - s'^q')]. (4.6) 
In the above equation Res[/(co')] stands for the residue of the function / at u. In the language 



of BPY [15|, the first two terms of the two spectral densities are the pole terms. The next 
two are the cut terms as they contain q and x; variables which will get integrated over 
before the k integration. There are, of course, other multiplicative factors which depend on 
k,x, and q which have not been expressly written down in Eq. ( [4.6|) . 

We now take the product of the two spectral densities as given in Eq. ([4.6[). Each 
combination of delta functions gives us a different cut of the diagram. For instance combining 
any of the pole terms from the two spectral densities gives us the cut shown in Fig. ^ (and 
thus, the resulting cross section for the second diagram ). This represents the process 
of gluon-gluon to e"^e~. This is a new process which to our knowledge has never been 
discussed before. The other possible cuts represent essentially finite density contributions 
to other known processes of dilepton production: two loop contributions to qq 7*, one 
loop contributions to qg 57*, and qq — > qq'j*. 

In the following we focus exclusively on the first process i.e., gg ^ 'y*. The reasons for 
this are twofold. The first is the fact that this process is unique and resembles nothing at 
zero density; all the other processes are extra contributions to processes which are known 
and have been calculated up to two- loop level at finite temperature and zero density |T^-|TB[. 
Secondly, all the other diagrams have at least one internal gluon line which is weighted by 
Bose Einstein statistics, and hence, these diagrams will display infrared divergence. To cure 
this divergence one can replace the bare gluon propagator with a resumed HTL propagator. 
These rather involved calculations deserve a separate treatment. In contrast, all the internal 
lines of the first process are fermions and hence show no infrared divergence. 
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V. THE CALCULATION 



We now apply the BPY formula with only the pole-pole terms from Eq. (|4.6| ). The 
various delta function combinations that are encountered can be schematically written as 



S{E — UJ — Lj')piP2 = S{E — u — u') 



5{iu - A:)Res[/i(cj = k)]5{iu' - A;)Res[/2l 



UJ 



+ 6{iu - k)Res[fi{iu = k)]6{iu' + k)Res[f2iuj' = -k)] 
+ 6{u; + k)Res[fi{uj = -k)]6{uj' - k)Res[f2{uj' = k)] 

+ d{uj + k)Res[fi{uj = -k)]6{uj' + k)Res[f2{uj' = -k)] 



(5.1) 



Note that only the first term (i.e., gluon-gluon annihilation) survives. In this term 
u = u' = k = E/2. Thus, this contribution to the discontinuity of the photon self-energy 
may be given as 

Ank'^dqdq'q'^q'^dede' sin 9 sin 6' 



Disc[n^] 



64 
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■'Q/37 



(27r)' 



(g - k)^^q^J]^{s2SiAn{Eg, gg/x) - S2S3An{Eq_k, Sip)} 
{k + {siEq_k - S3Eg)){k - {siEg_k - Sa^J) 

il' - f^)s[,a'l's'^,p'Q's'„YW2s'iAn{Eg>, s'^fi) - s^s'3An(E,/_fc, s»} 
{k + {s[Eg,^k - s'3Eg,)){k - {s[Eg,^k - s'2E,,)) 



(5.2) 



As stated before, the multiplicative factor C'^p^'^ is the result of the contraction of Lorentz 
indices appearing in the two traces of 6 7 matrices and the intervening factors of the metric. 
Essentially, this serves the purpose of the average over initial spins. We basically now have 
four terms to integrate (there is no mixed term between (g, 6') and {q',0') ). Each term is 
classified by its fermionic distribution function. Now we may shift the variable q ^ q + k 
in the second and fourth fermionic distribution functions, and the 6 and 6' integrals become 
simple and can be done analytically. The two q and q' integrations are then done numerically. 

The differential production rate for pairs of massless leptons with total energy E and 
and total momentum p = is given in terms of the discontinuity in the photon self-energy 
as 



IS 



dW 
lEdfp 



(p = 0) 



a 



1 



1 



1 



127r3 ^2 1 - e^/^ 2 



Discn;;(o). 



(5.3) 



Where a is the electromagnetic coupling constant. The rate of production of a hard lepton 
pair with total momentum p = at one-loop order in the photon self-energy (i.e., the Born 
term ) is given as 



dW 
lEdFp 



(p=0) 



a 

6^ 



h{E/2- p)h{E/2 + p). 



(5.4) 



As mentioned before, initial temperatures of the plasma formed at RHIC and LHC have 
been predicted to lie in the range from 300-800 MeV [ll3|,p^ . For this exploratory calculation 



10 



we use a conservative estimate of T = 400 MeV. To evaluate the effect of a finite chemical 
potential we perform the calculation with two extreme values of chemical potential fi = O.IT 
(Fig. H) and /i = 0.5T(Fig. ^ ^^j. The calculation, as stated before, is performed for three 
massless flavours of quarks. In this case the strong coupling constant is (see ||18|) 

"^^^^ = 27ln{T/50MeVy ^^'^^ 

The differential rate for the production of dileptons with an invariant mass from 0.5 to 
2.5 GeV is presented. On purpose, we avoid regions where the gluons become very soft. In 
the figures, the dashed line is the rate from tree level qq (Eq. ( ^.4|) ); the solid line is that 
from the process gg e~'"e~. We note that in both cases the gluon-gluon process dominates 
at low energy and dies out at higher energy leaving the qq process dominant at higher energy. 



10" 



10 



10" 



10" 



10 




FIG. 4. The differential production rate of back to back dileptons from two processes. Invariant 
mass rmis from 0.5 GeV to 2.5 GeV. The dashed line represents the contribution from the process 
qq — > e+e~. The solid Hne corresponds to the process gg — > e+e". Temperature is 400 MeV. Quark 
chemical potential is O.IT. 



VI. DISCUSSIONS AND CONCLUSION 



In this paper we have performed a calculation to estimate the effects of a non-zero quark 
chemical potential on the intermediate mass dilepton spectra. We have found that a new set 
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0.5 1.0 1.5 2.0 2.5 

E (GeV) 

FIG. 5. Same as Fig. | but with fi = O.ST 



of diagrams may become important at finite temperature and finite density. These diagrams 
lead to a new contribution to the photon self-energy at the 3-loop order. There are various 
cuts to this diagram. Most of these result in finite density contributions, and/or higher 
order excess contributions to well known processes. One of the cuts however represents 
an entirely new process. The contribution of this diagram to the differential production 
rate of back-to-back (or p = 0) dileptons is estimated. This is then compared with the 
contribution emanating from the tree level process of qq annihilation. The rate from this 
new process is found to be larger than the simple tree level rate by orders of magnitude 
between a dilepton invariant mass of 0.5 to 1.5 GeV. The reasons for this large magnitude 
are many. The gluon-gluon diagram is enhanced by the Bose-Einstein distribution function 
of the gluons. Also, at lower energy the An factors which represent the difference between 
the quark and antiquark distribution functions are comparable in size to the distribution 
functions themselves. There is, also, enhancement from the larger color factors of the gluons. 
One also notes that as the energy (invariant mass) is increased, the difference between the 
two distribution functions (An) decreases rapidly, i.e., as higher energies the distributions 
of quarks and antiquarks becomes less sensitive to the chemical potential. This is the main 
reason behind the sharp drop of the differential rate compared to the differential rate from 
quark- antiquark annihilation. 

These results demonstrate the importance of these finite fi processes on the dilepton 
spectra emanating from a quark gluon plasma with a quark-antiquark asymmetry. It is 
simple to note that this diagram is most sensitive to gluon number. The early stages of the 



plasma have been predicted to be gluon dominated pT| . The contribution from this diagram 
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should clearly shine in such an environment. 

The treatment in this rather exploratory work should, and will, be improved upon. Our 
goal here was simply to establish the existence of a signal. One may have desired, for 
example, that the calculation be extended to arbitrary p. However, similar extensions, even 



in the HTL approximation, are known to be rather involved |T2[. There are also the other 
diagrams which have yet to be computed. These however suffer from the defect of having an 
internal boson line which will show an infrared divergence. This bare propagator will have 
to replaced by a resumed HTL propagator, in the event that the momentum flowing through 
it becomes very small. These aspects, along with others, will be addressed elsewhere. 
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VIII. APPENDIX A 

Notation. 

Our notation is categorized by the explicit presence of an apparent Minkowski time 
= —IT and a momentum g° = i{2n + 1)tiT + /x. Our metric is (1, —1, —1, —1). For the 

case of zero chemical potential our bosonic propagators have the same appearance as at zero 

temperature, i.e.. 

The Feynman rules are also the same as at zero temperature, with the understanding that 
we replace the zeroth component of the momentum by Eq. ( p^.2| ) for a fermion and by an 
even frequency in the case of a boson. One may, in the case of zero chemical potential, relate 
this to the familiar case of reference |]T0[ by noting that 



1 -1 

^(^) = T^2 ^ = 7 I I 12 = -^E{^n, q), (8.2) 

where A£;(c<j„,g) is the familiar Euclidean propagator presented in the literature ( pO| , 
||TT|). One may immediately surmise the form of the non-covariant propagator A(|g|,a;'^), 
the Fourier transform of which is the covariant propagator. 
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—I 



-1(3 



X 



0^ 



(8.3) 



In the presence of a finite chemical potential our full fermionic propagators become 

J 

where, 



•4) 



(8.5) 



In the above two equations g° = i{2n + 1)ttT + fi. 



IX. APPENDIX. B 



Derivation of Eq. |3.4| . 

As stated in section III we start by substituting Eqs. ( |XT| , |H^ , [H7^ ) in Eq. (|2.1j). Substi- 
tuting in the first diagram we get, 



{q + p-k)s,aq-s2p{(l + P) 



S37 



^-iS2{Eg+S2ll)x^^-iSl{Eq^p_i:+Slll)x1^-iS3{Eg+p + S3ll)x'^ 



-03 



dec dec 2 dx 2 



(9.1) 



In our formalism x° = —ir. Also, in Eq. (|2.1|) the zeroth components of the momenta are 
replaced by time derivatives e.g.. 



1o = '-no+ 



with the derivative acting on the e~*'^'^°~'^)^2. We then perform a by parts integration w.r.t. 

after which the derivative acts on the e-«'52 (£^9+52^)2:2 term. After some algebraic ma- 
nipulation and using the definition that 5^20 = S2Eg, we obtain Eq. (|9.1|). We have the 
identity 



1 ^ g.,0(,o_,o+,0) ^ ^^^^ _ ^ ^ _ + ^3 _ 



(9.2) 
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Using the above identity we perform the Matsubara sum and then evaluate one of the r 
integrals using the two delta functions. Noting that = 1 we get 



27r3 2 



■Tr [67] J2 



S37 



Sl,S2,S3 



-if3 



-ix'^{S2Eq + S3Eg+p) -ix1{siEq^p_^.-S3Eq+p) 



0(r2 - Ti) - 0(ri - T2)e 



l3{-s-3Eg+p+n) 



(9.3) 



Noting that e^(-'*3^9+f+'') =- - 

fs3{Eq+p,S:jll) 

two Heaviside theta functions to obtain, 



/ q+p, 3fj-) niay integrate X2 making explicit use of the 



rp^pu 



'ah 



27r3 2 

,0^1,0. 



Tr [67] 



S37 



Sl,S2,S3 



-il3 



^Eqj^p_.}^EqEqj^p 
^-^4i^^+s,Eq,p.,-S3E,)J^^^j^^^^_^^ S^^)fjEq, S2f^) 

^i{-il3)ljP -S2Eq-S3Eq+p) _ ^i{x^){p° -S2Eq-S3Eq+p) 



dxl 



+ f_s^{Eq+p, Ssfl) 



i{p° - S2Eq - SsEq+p] 

^i{x1){p'^ -S2Eq-S3Eq+p) _ ^' 

i{p° - S2Eq - SsEq+p) _ 



(9.4) 



We note that, 



^-f3(s2Eq+S3Eq+p) _ ^-t3s2{Eq + S2fi) ^-t3s3{Eq+p-S3p) 



^ f-SijEq, S2fi) f'_,.^{Eq+p, S3/i) 
f^^{Eq,S2fl) fs,{Eq+p, Ssfi) 

Substituting the above equation in Eq. ( |9.4|) , and noting that fs — fL 
the integration to obtain 



sal 



(9.5) 

s, we may perform 
1 



2^' 2 ''^^^.Xs, ^Eq+p-kEqEq^p pO - S2Eq - S^Eq+p 

fsAEq+p-k-, Slfi)f'_s.^{Eq+p, g3/i) — f'_^_^{Eq+p_k, Sifl) f+sajEq+p, ^3/1) 

+ SiEq^p^k — S^Eq^p 
fsAEq+p-k, Slfl)fs2{Eq, S2H) - f'_s,{Eq+p-k, Sl^l) f'_^^{Eq, gg/^) 



-S2 



JsiK^q+p-k, aifJ,)Js2KJ^q, -^2/^; - J _s^K^q+p-k, a l^^ ) J _ s^K-^^ q ^2^') /q ^^ 

' kO-p^ + S,Eq+p_k + S2Eq J" 

Expanding the /'s in terms of Fermi-Dirac distribution functions, we finally obtain the full 
expression for T^^'^ as (note: S3 has been changed to —S3 ) 

^ - ^ ^ 



-S37 



1 



27r3 2 



Sl,S2,S3 



'■^q+p-kEqEq^p — S2-Eg + S-^Eq^p 

{Sl + S3)/2 - S2.n{Eq+p_k - Sifi) - Sin{Eq+p - S3/i) 



53 



(si + S2)/2 - S2n{Eq+p_k - sin) - Sin{Eq + Sa/i) 

k^ -p° + SiEq+p^k + S2-Eg 



-S2 



(9.7) 
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We, now, perform the same set of manipulations for the other Feynman diagram of Fig. 
to obtain 



T 



ab 



27r3 2 



Tr [67] J2 



-S37 



1 



- S3 



Si + S3)/2 - S3n{Eg^p+k - sifi) - sih{Eq_p - s^fi) 

+ SlEq_pj^k + S3^g-p 

[si + S2)/2 - S2n{Eq_p+k - sifi) - Sin{Eq_p + S2^) 



S3E, 



q-p 



k'^-p^ + SlEg^p + k + S2Eg 



(9.8) 



In the above two equations we have Tr [67] = Tr [7'^7'^7^7''7'^7"] in the first equation and 
= Tr [7'^7''7^7^7'^7"] in the second equation, note that both these factors are the same. The 
only difference between the two equations is the sign of q. We now set q — > —q in the three 
integration of the second term. As a result {q — p + k)sj^a —{q + p — k)-sia,q-s2fB ~^ ~1s20 
and (g— p)_s37 — > — (g+p)<j3^. This is then followed by setting si — si, S2 —S2, S3 —S3 
in the summation over Si, 52,^3. With these changes to Eq. (|9.8|) , we may now add Eqs. 
(|9.7D and (|9.8|), and finally change (s2, S3) — > (— S2, —S3) to give Eq. 



X. APPENDIX. C 



Interdependence of the T^p" integrals 

In this appendix we will analyse the (potentially 64) terms generated by 



J 



0/37 



Tr [7A'^/3^p^'5^i'^" j 



:io.ii 



We note that the Lorentz indices affect, specifically, only part of the integrand of J'^l^i ^ i.e., 
we may write 



= j ded(f) sine 



r^\9,<P)i dqq'j{q,9,ft) 



;io.2) 



As noted in section III, performing the integration sets 44 of the 64 terms to 
zero, the only surviving terms are those that carry (a,/3, 7) in some permutation of 
the combinations (0,0,0), (3, 3, 3), (0,0, 3), (0,3, 3), (0 integration gives 27r), and those of 
(0,1,1), (0,2, 2), (3, 1,1), (3, 2, 2), (0 integration gives vr). The /i dependence in J only distin- 
guishes between terms with /i = and /i 7^ 0. On performing the integration we get 



J 



a/37 



dd sin d 



/"^"W / dqq'j{q,9,fi)]F{a,(3,j) 



;io.3) 



where J{q, 6, /x) = Jof^/^.o + Ji(l - 5;,,o), and 
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F{a, P, 7) ^n[{ga,igp,l + ga,2gp,2){g'y,0 - ^7,3) + {9a,l9j,l + 9a,2gj,2){gp,0 - 9(3,3) 

+ {30,l9^,l + 9/3,29y,2){9a,0 " 9a,3)] + '^7^ [{9 a, 9/3,0 + 9a,39/3,3){9-/,0 " 9'r,3) 
~ 9a,09l,09P,3 + 9a,39l,39(3,0 — 9l3,09t,09a,3 + 9/3,39f,39a,0 ■ 

We now look at the integrands of the 9 integration: note for example that 
/ooo = S1S2S3; /°°^ = siS2Cos^; = sicosOsa. 
The reader may easily verify that this implies that 

j-003 



^030 
S2 



Also note that 



which implies that 



r033 



Si cos^ 9; I 



oil 



r022 



Si sin^ 9, 



oil 



J 



022 



J033 ^ 2 J^^^ = 



Oil 



J' 



000 



S2S3 



Following the above method we may also prove that 



J 
J 



303 


— — ^330 






S3 




311 


_ j322 




131 


_ j-232 _ 


^113^ 


333 


+ 2 J^i^ = 


^300 






S2S3 ' 



J 



223 



:io.4) 



;i0.5) 



:io.6) 



The above 12 conditions reduce the number of independent J'^^^^s to 8. Thus, one only 
has to evaluate these 8 terms and the others may be evaluated by the above mentioned 
conditions. 



XI. APPENDIX. D 

Derivation of the BPY formula. 

In this section we present a discussion on the BPY formula. The basic aim is to evaluate 
the quantity 



Disc 
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where p^,k^ are both discrete frequencies. As stated in |TT|] the sum over the Matsubara 
frequencies can be converted into two contour integrals over a complex ( Eq. (3.39) of 
m ), i.e., 



1 



,1 



1 



too— e 
ioo+e 



[11.1) 



Now fi{z) and f2{z) both have discontinuities on the real axis, which consists mostly 
of residues and branch cuts. Let the sum total of all such discontinuities be set equal to a 
spectral density p{uj), i.e., 

2nip{uj) = J25{uj-pi)27riRes[f{z = to)] +J2e{uj - I Q {u j - u)Disc[ f{z = u)], (11.2) 



where pi represent the location of the poles of f{z) on the real axis ,and uj and Ij represent 
the upper and lower bounds of the branch cuts on the real axis. The presence of p° in 



the arguments of only one of the functions in Eq. (|1 1 . 1| ) separates the discontinuities of 
the two functions. We have tacitly assumed that both these functions contain integrable 
discontinuities. Thus the contour integration can be reduced to a sum over residues and 
discontinuities of the two functions separately. Let us look at the first integral. The function 
fi{k^) {k^ now complex ) has discontinuities at = u;( real). The spectral density pi{uj) is 
defined as 

2iTipi{uj) = [fi{z = UJ + it) - fi{z = u - it)], 
thus the total contribution of this term to the clockwise integral is 



dupi{uj)f2{p^ - uj) 



1 1 

2 ^ 1^- 



The function g(kP) = f2{p^ — k^) has discontinuities at k^ = p^ 
density now, is 



uj{ real). The spectral 



g{z = p — u + it) — g{z = p 



UJ — tt] 



= f2{z = uj-it)-f2{z = uj + it) = -2mp2{uj). 
Thus the contribution from this term to the first clockwise integral is 



duj'p2{uj')h{p'' -Uj') 



+ 



2 ■ Qfi(p'>-i^') _ 1 



Performing the above mentioned operation for the other contour integral we get. 



S 



dujpi{uj)f2{p^ - uj) 

/ dwpi{uj)f2{p^ - ^ 
Jo 



- 1 



+ 1 duj'Mp''-uj')p2{uj') 
'o 



+ 



1 



1 1 

2 ^ e^'^ - 1 



+ / duj'Mp'-Uj')p2{uj') 



1 

2 ^(p"-'^') - 1 



2 ■ e^(p"-^') - 1 
1 



11.3) 
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The two lines in Eq. ( [11. 3| ) display the contributions from the two separate contour integrals 
of Eq. ( 11. 1| ). These may be combined together to give, 



+ 



dujpi{uj)f2{p^ - uj) 

oo 



1 1 

2 ^ IP^- 



1 



1 1 

2 + e/3(p°-^') 



(11.4) 



Now, we start the operation of analytically continuing ]P onto the real axis. Note that 
the Bose factor of the second term will have a discontinuity when ]P is analytically continued 
onto the real axis. Thus the correct analytic continuation is achieved by dropping the p° 
from the expression (as exp/^p*^ = 1 when p*^ = —i2nnT). We observe that in the first 
term of Eq. ( [11. 3| ) the only factor that may have a discontinuity is /2(p° — uj) and this 
only happens when — uo = uj'{a. real number). This is only possible when E ± ie. 

Taking the discontinuity across the real axis in both terms of Eq. ( [11.3|) , and combining 
them together we get. 



Disc[5(E)] = -2m 



duj 



duj' pi{uj)p2{uj')5{E — UJ — uj')< 



;ii.5) 



We observe that the term in the curly brackets may be simplified: 

1 1 ef^-l 



(e/3- - l)(e-/3-' - 1)^ 



- 1 e-/5'^' - 1 
thus we get back the BPY formula, i.e., 

DiscTj] /i(A;°)/2(/ - A:°) = 2m{l - e^/^) / du / du' 

n{ijj)n{uj')5{E — uj — uj')pi{uj)p2{uj'). 



ko 



(11.6) 
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